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ABSTRACT 

Aims. Bipolarity in proto-planetary and planetary nebulae is associated with events occurring in or around their cores. Past infrared 
observations have revealed the presence of dusty structures around the cores, many in the form of disks. Characterising those dusty 
disks provides invaluable constraints on the physical processes that govern the final mass expulsion of intermediate-mass stars. We 
focus this study on the famous M2-9 bipolar nebula, where the moving lighthouse beam pattern indicates the presence of a wide 
binary. The compact and dense dusty core in the center of the nebula can be studied by means of optical interferometry. 
Methods. M2-9 was observed with VLTI/MIDI at 39-47 m baselines with the UT2-UT3 and UT3-UT4 baseline configurations. These 
observations are interpreted using a dust radiative transfer Monte Carlo code. 

Results. A disk-like structure is detected perpendicular to the lobes and a good fit is found with a stratified disk model composed 
of amorphous silicates. The disk is compact, 25x35 mas at 8yum, and 37x46 mas at 13//m. For the adopted distance of 1.2 kpc, the 
inner rim of the disk is ~15 AU. The mass represents a few percent of the mass found in the lobes. The compactness of the disk 
puts strong constraints on the binary content of the system, given an estimated orbital period 90-120yr. We derive masses of the 
binary components between 0.6-1.0Mq for a white dwarf and 0.6-1.4Mq for an evolved star. We present different scenarios on the 
geometric structure of the disk accounting for the interactions of the binary system, which includes an accretion disk as well. 

Key words, planetary nebulae - post-AGB stars - infrared interferometry - disks 



1. Introduction 



Intermediate mass stars (I-8M0) undergo extreme mass-loss during the late stages of their evolution. The mechanisms that dom- 
inate those stellar outflows are also responsible for the shaping of the ejecta. Hints of asymmetries are numerous in many stars 
characterizing these late stages. Almost all of planetary and proto-planetary nebulae (Sahai et al. 2007) observed by the Hubble 
Space Telescop^jexhibit asymmetries in their shape, while bipolarity and equatorial dusty disks are very common as well. 

Most proposed shaping mechanisms for bipolar nebulae involve binarity (Zijlstra 2007 ), since a source of angular momentum is 
required to drive the ejecta in directions perpendicular to the orbital plane. Interaction between the ejecta and the orbiting companion 
can lead to the formation of torii, disks and jets, the latter perpendicular to the binary system's plane. During common envelope 
evolution, accretion disks can be formed leading to the ejection of jets, thus forming bipolar nebulae (Soker & Livio 1994). Disks 
can be formed during wind accretion preceding common envelope stage as well (Mastrodemos & Morris 1998 ), but are too low 



mass to be the direct cause of the shaping of the older nebula ( |Chesneau et al. 2007 [ ). The disks appear to be a remnant of the shaping 
process. 

Disks are rather small in size (a few milliarcseconds), compared to the surrounding nebulae (many arcsecs), and their dust re- 
emits stellar radiation in the infrared. Observations of such small structures near evolved stars can only be performed with the use of 
high-angular resolution techniques like the infrared Very Large Telescope Interferometer (VLTI), i.e. planetary nebulae: Chesneau 
et^p006||2007l>;|Lagadec et al.|l|200o); VLTP star : jChesneau et al.|([2009); Mira stars: |Tatebe et al.| ( p006l [2007| |2008)l; Carbon " 
stars: |Ohnaka& Boboltz (2008); Ohna ka et d.|l|2008);|Sacuto et al.|(|2008|);|Deroo et al.|(|2007b) ; OH/IR stars: |Deroo et al.|(|2007a|b); 
Chesneau et al.|(|2005b[); po st-AGB stars: Deroo et al. (2006i; Matsuura et al. ( 2006| l; proto-planetary nebulae: Murakawa et al. 
( |2008| l; |Ohnaka et al.| ( |2006) . 



1.1. Twin-jet nebula M2-9 

Two of the best-studied examples of elongated and axisymmetric bipolar planetary nebulae are M2-9 and Menzel 3. They present 
tightly pinched waists and symmetric lobes. They are also spectroscopic twins at visual and IR wavelengths ( |Smith|2003" Smith & 



Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere, Chile, ESO N: 079.D- 
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http://www.astro.washington.edu/users/balick/WFPC2/index.html 
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Table 1. Observing log 



OB 


Time 


Base 


Projected baseline 








Length 


PA 








[m] 


[degrees] 


M2-9_l 


2007-04-11 


U2-U3 


40.2 


107.3 


M2-9_2 


2007-06-28 


U2-U3 


47.3 


126.0 


M2-9_3 


2007-08-28 


U2-U3 


39.7 


44.0 


M2-9_4 


2008-03-28 


U3 -U4 


45.1 


39.0 



Calibrators: HD 116870 K5III 2.44 ± 0.12mas, HD 163917 G9III 2.75±0.11mas, HD 167618 M3.5III 11.33±0.1mas, HD175775 G9III 

3.26±0.23mas, HD 152820 K5III 2.57±0.34mas 



Gehrz 2005 ). This paper is focused on M2-9. The lobes of the nebula extend up to 22" at each side of the bright central core. Schwarz 



et al. (19975 found that the nebula has polar knots extended up to 52" from the central source (core). The core is dominated by Ha 



emission scattered by dust in the visual, as well as free-free emission in radio (diameter < 3"). 

A lighthouse beam is seen inside the bipolar lobes and a series of observations at different epochs revealed that it is a by-product 
of the mechanism responsible for the shaping of the planetary nebula. The lighthouse effect was witnessed by knots in the nebula 
known as early as 1952 (All en & Swings|1972| l, was monitored from 1989 to 2007, revealing a period of ~100 years. This feature 
can be created by orbital motion and is the best evidence that the core engulfs a binary system. In addition to that, the high-excitation 
lines seen in M2-9 suggest that the illuminating source is a hot and compact star, i.e. a white dwarf (WD) ( Schwarz et al.| 1997 Lim 



& Kwok 2000; |Livio & Soker 2001 ), not to be confused with the primary star whose ejecta now constitute the nebula. In such a 



scenario, dust originating from the primary will settle in the binary system's orbital plane in the form of a disk, torus or spiral. 
The distance of M2-9 is uncertain. Many authors quote a distance of 1,000 pc in the absence of a tight constraint (|Hora & 



Latter|[T994l |Livio & Soker||2001| |Smith & Gehrz |[2003] l. |Schwarz eTIiLl (fl"997) propose a much closer distance of 640±100pc, 



based on expansion parallax measurements of the outer lobes. We derive a larger distance in this paper, confirmed by R. Corradi 
(private communication) who reanalyzed the data of |Schwarz et al. ( 1997| l. For all the mass estimates described below, the assumed 
dis tancapjwas diffe rent from the one presented in this paper. 



Bachiller et al. 



Mq for the entire nebula, but |Lenzuni et al. ( 1989) estimated a higher total 



1990) mount of molecular gas in the nebula h igher than 10 3 Mrr>. 

gaseous mass of 0.21 Mq and mass of dust 2.1 x 10" 
dust mass by IRAS observations (~ 5.4 x 10~ 3 M©). 

The core of M2-9 was resolved using millimetric interferometry (Plateau de Bure, 5"x3" beam), showing that the molecular 
gas is distributed in a ring-like structure extended about 6" E-W. A kinematical age ~2,100 years was derived from the expansion 
of the CO gas in the nebula ( Zweigle et al.|1997] >. 

Observations by Smith & Gehrz ( 2005 \ within the last decade showed that in the core the mass of warm gas is 9.2 x 10~ 4 Mq and 
that of warm dust 4 X 10~° Mq, while in the lobes the mass of cold gas is 0.78 Mq and that of cool dust is 3.4 x 10~ 3 Mq. 

We present observations of the core of M2-9 in the near- and mid-infrared (Section 2). We have detected a dusty disk in the 
core with the use of high-angular resolution techniques. The geometric constraints of the disk were defined with the use of radiative 
transfer models. Our results are presented in Sections 3 & 4. A detailed discussion can be found in Sections 5-6, where we elaborate 
scenarios on the disk's structure. In Section 7, we compare our findings with the disk of Menzel 3. Our conclusions are in the final 
section. 



2. Observations 

2.1. VLTI 

M2-9 was observed in the N band (7.5-13.5/zm) with MIDI ( |Leinert et al.|2003|[2004| l, the mid-infrared recom biner of the VLTI 
with the use of only two Unit Telescopes (UTs). We used a typical MIDI observing sequence as described in Przyg odda et al.| 
(2003). MIDI can make single-dish acquisition images with a field-of-view of about 3" with a spatial resolution of about 0.25" at 
8.7//m, and provides a flux calibrated spectrum at low spectral resolution (in this case R=30) and several visibility spectra from the 
sources (Chesneau et al. 2005a b), i.e. spectrally dispersed information on the spatial extension of the source. Only four visibility 



measurements could be obtained, due to bad weather conditions. 

The observations of M2-9 with MIDI were performed in April, June and August 2007 and March 2008 in the so-called 
SCIJ-'HOT mode, meaning that the photometry of each telescope is recorded simultaneously to the fringes. The errors, includ- 
ing the internal ones and the ones from the calibrator diameter uncertainty, range from 5% to 10%. The accuracy of the absolute 
flux calibration is about 10-13%. For the flux calibration, HD 163917 (F^m = 16.54 Jy) was used. A Cohen template for a K0III 
star was used (HD 180711, F^m = 21.35 Jy). The log of the observation is given in Table [T] We used two different MIDI data 
reduction packages: MIA developed at the Max-Planck-Institut fur Astronomie and EWS developed at the Leiden Observatorjj^] 



All mass estimates with an assumed distance of lkpc, except for|Lenzuni et al. i 1989} (2.37kpc) 
Available at http://www.strw.leidenuniv.nl/~nevec/MIDI/index.html 
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Arcsecond 

Fig. 1. Composition of HST/STIS (black) and NACO/Ks (red) images. North is up and East is on the left. The infrared emission 
contours follow the lines of the HST scattering in the visible, not only in the core but in a N-W region as well. The unresolved 
core is less than 0.1" in Ks. Contours from the center to low levels are 95, 10, 2, 0.8, 0.4, 0.2, 0.1, 0.05, 0.025 and 0.015% of 
the maximum, respectively. Noise appears in the last two lower levels and those reveal the extended nebula as it is evident on the 
northern and southern part of the image. The diffuse light is increasing smoothly from the 0.05 to the 10% level. 



2.2. NACO 



We observed M2-9 with the near-infrared adaptive optics instrument NACO fLenzen et al.|2 003 ; Rou sset et al.|2 003 ) as well, using 



three broad-band filters centred respectively at 2.18/im (Ks), 3.8yum (L') and 4.78//m (M')- The Ks data were taken during the night 
of 3 July 2006, and the L'-band and M'-band data on 18 August 2006. These observations were complemented with close-by PSF 
calibrations using HD155078 (F5V) in V and HD156971 (FUJI) in M'. We used the S27 (Ks) and L27 (L', M') camera mode to 
obtain a field of view of 28"x28"and the pixel scale of 27.1 mas per pixel. Jittering was used to remove the sky and for the M' 
observations chopping was also used. The reduction procedure followed for MIDI and NACO data of M2-9, is the same as the one 
described for Mz3 in Chesneau et al. (2007}. 

Our Ks, L' and M' observations of Mz3 (Chesne au et al.|2007| > and M2-9 with NACO revealed an unresolved core at diameters 
less than 0. 1" . Figure [T]is a composite image of a Hubble Space Telescope (HST) STIS/CCD image in the visible and our Ks-band 
image, the latter being an improved near-infrared view of the source compared to Aspin et al. (1988). One can notice a broad 



near-infrared emission around the core of M2-9 (diameter< 5") and some enhanced emission in the lobes as well (~10" from the 
core; more prominent at N-W). These are attributed to scattering. Near-IR photometry gave m/,<=4.12 and m^=2.5 and these are 
included in Fig. [6] 



2.3. ISO 

We have reduced and analysed archivecQdata from ISOCAM ( Cesarsky et al. 1996 1, the Infrared Space Observatory's (ISO) infrared 
camera. The satellite was mounted with a 60cm Ritchey-Chretien Cassegrain telescope. Ten exposures have been recorded in 1996 
for a range of wavelengths from 3.0 to 14.9/im with both short-wave (S W) and long- wave (LW) filters (bandwidths stated in Table[2]i. 
Each frame is 32x32 pixels with an effective field-of-view 45"x45" and the spatial scale of each pixel is 1.5". The standard software 
CIA was used to reduce and analyse the data (Table[2} in the procedure described in the manual. Photometric errors were consistent 
for each detector. In ISOCAM the physical detector pixel size is 1" and the Airy disk pattern increases with wavelength from 1" to 



4 http://iso.esac.eso.esa.int/ida/ 
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Table 2. ISOCAM photometry 



[j"m] 


BW [/an] 


detector 


Flux [Jy] 


o-flux [mJy] 


3.0 


1.00 


SW3 


4.6 


9.5 


3.7 


0.55 


SW6 


7.9 


9.6 


4.5 


0.30 


SW7 


13.3 


9.5 


4.5 


1.00 


LW1 


16.1 


16.7 


6.0 


1.00 


LW4 


27.6 


16.7 


6.8 


0.50 


LW5 


32.1 


16.4 


7.7 


1.50 


LW6 


37.2 


16.6 


9.6 


2.20 


LW7 


34.7 


16.7 


11.3 


1.30 


LW8 


42.3 


16.7 


14.9 


2.00 


LW9 


50.6 


16.6 



3.5" (see[Blommaert et al. 2003 p. 44) thus the core of M2-9 was not resolved by this instrument. The broad infrared emission from 
the core did not permit the detection of any other significant structure in its vicinity either. 

An infrared spectrum has also been retrieved from ISO Short Wavelength Spectrometer (S WS) using short- and long-wavelength 
1996] >. The aperture area varies with wavelength range and grating: it is 14"x20" for 2.38-12.0/im, 
for 12.0-27.5/zm, 20"x27" for 27.5-29.0//m and 20"x33" for 29.0-45.2^01 (see FigM. 



gratings (Kessler et al. 
14"x27' 



2.4. Spitzer 



MP Et i HSI [hSMHJ 




Fig. 2. Positioning of the Spitzer/IRS slits over a F656N exposure of M2-9 (Hubble Space Telescope Legacy Archive). Both slits 
were centered on top of the northern lobe, the Short-High (green; 4.7"xll.3") at a PA= 13 1.2° and the Long-High (red; 1 l.l"x22.3") 
at a PA=46.1°. North is up, East is left. The underlying Ha image comes from 1996. By the time that Spitzer observed M2-9, the 
underlying structure of M2-9 was rather different, and the nebular structure inside the green aperture had changed. 



Calibrated data were extracted from the Spitzer Space Telescope^] Archive. M2-9 has been observed with both the InfraRed 
Array Camera (IRAC) and the InfraRed Spectrograph (IRS). Unfortunately, the core's strong infrared emission has saturated IRAC 
rendering those observations useless. 



5 http ://s sc . spitzer.caltech .edu/ 
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High-resolution spectra (Fig. [3} by IRS revealed important information on the nebular dusty and gaseous chemistry. Two slits 
were positioned over the northern lobe and avoided the core (Pig. |2j: Short-High slit (9.9-19.6yum) is 4.7"xl 1.3" at PA=131.2° and 
Long-High slit (18.7-37.2/xm) is 11. 1"x22.3" at PA=46.1°. 

2.5. AKARI 



We include mid-infrared data from t he AKARI satellitepj lRC Point Source Catalogue archive ( Murak ami et al.|2007| l. M2-9 has been 



observed with the InfraRed Camera ( [Onaka et al. 2007) > at two infrared bands, namely S9W (6.7-1 1.6jum) andL18W (13.9-25.6/mi) 



with respective effective wavelengths at 9 and 18 /mi (see Fig. |6j 

3. Results 

3.1. Spectral Energy Distribution 

The integrated flux, as detected by MIDI (field-of-view: ~0.4" with UTs), in Fig. [6] clearly shows a 20% difference between the 
MIDI and the ISO spectra. Since the rectangular aperture of ISO (14"x22") was positioned over the central area of M2-9 (excluding 
the largest part of the lobes) and it increased in size for bands larger than 12yum (including part of the lobes near the central source), 
it is expected that ISO detected a more extended source than MIDI. The fact that both spectra have similar shape and only slightly 
different flux levels, indicates that they both detected almost the same compact source. 

Both ISO and MIDI spectra, shown in Fig. [6] have been dereddened using the Savage & Mathis ( 1979) law for Ay — 2.5 ( |Torres- 
|Peimbert & Arrieta|1998| l. 

The two Spitzer slits were positioned on top of the northern lobe, but they have probed different parts of the nebula. In more 
detail (Fig. [3): 

- Short-High: 

The spectrum is very weak compared to the continuum of both MIDI and ISO, which were positioned over the core, but 
there are many interesting features. There is a broad PAH emission at 11/mi and much fainter broad regions at 12.3 and 
12.8//m. [Nell], [Nelll], [Fell] and [SIII] lines are seen at 12.8, 15.6, 17.9 and 18.7/im, respectively. No other PAHs have been 
detected near the core by ISO, whereas it is unclear whether the 11.4/im emission at 10% of the continuum of MIDI spectrum 
originates from PAH (Fig. [4]). This suggests that the 1 lyum feature resides in the lobe. 

- Long-High: 

There is an indication of a broad feature at 33.5//m (forsterite) but the neighboring broad shallow absorption from 30-34//m 
suggests that they may both not be real. Prominent gaseous emission lines are also present: [Felll] (22.9/mi), [OIV] (25.9/mi), 
[SIII] (33.48/mi) and [Sill] (34.8/im). All those require a high excitation state and would be expected near a white dwarf. 

3.2. Interferometry 

Our VLTI observations of the core of M2-9 reveal the presence of a dusty flat-like structure. Visibilities decrease as the PA. of the 
projected baselines increases (Pig. [5j, suggesting an increase of the dust opacity. Thus, MIDI has probed the existence of a flattened 
structure along the equatorial plane of the nebula. 

The core is significantly resolved by the 40-50m baselines, with visibilities of the order of 0.1-0.2. Our observations showed 
that the flattened structure is compact, typically 25x35 mas at 8//m, and 37x46 mas at 13//m (Fig. [5}. The dusty structure is more 
elongated along the planetary nebula's equatorial plane, compared to the direction of the lobes. Therefore, the disk's geometrical 
size can be constrained (see Sections 4 & 6). 

Subtracting the continuum from each spectrally dispersed visibility, one may enhance subtle visibility fluctuations that could not 
be distinguished otherwise. Significant signatures of forsterite (Mg2SiO/0 are observed as small (~ 5 — 10%) dips in the visibility 
curves (Fig. [4]) at 9.8 and 11.4/mi ( Molster et al.|2002 i. We confirmed the validity of this, since both features were present in the 



uncalibrated visibilities and they were absent in the calibrator data. In addition, the peaks nearly match the emission curves of 



forsterite grains measured in the lab by Koike et aL] (2000 2006). The visibility decrease at typically 10% shows that the forming 
regions of these features are resolved and thus slightly more extended than the continuum and amorphous silicate regions (as already 
observed in Young Stellar Objects, van Boekel et al.| ( [2004| l). 



4. Physical parameters of the disk 

In this section we present the modelled disk and the constraints introduced by the fitting procedure. 
4.1. Model 

The 3D radiative transfer code MC3D, which is based on the Monte Carlo method and solves the radiative transfer problem self- 



consistently, has been used in this work ( Wolf et al. 1999 1 Wolf 2003). It simulates the temperatures produced by dusty configurations 



6 http://www.ir.isas.jaxa.jp/AKARI/ 



6 



F. Lykou et aL: A disk inside the bipolar planetary nebula M2-9 



and creates as observables spectral energy distributions and wavelength-dependent images of the dusty environments, as well as 
polarisation maps. The density law used in this model is 



P(r,z)=p (^) exp -i(^) 



(1) 



h(r) = ho\ — 



(2) 



where r is the radial distance in the disk midplane, R* is the stellar radius, j3 is the vertical density parameter, a is the density 
parameter in the midplane and ho is the scale height at a given distance from the star (Shakura & Syunyaev 1973| Wood et aLp 002 ) 
. We have used the distribution of Mathis et al. ( 1977 1 for the grain sizes: n(b) oc b" 3 ' 3 , where b is the grain size, assuming that the 
dust grains have spherical symmetry. 

The dust residing around evolved stars may be oxygen- or carbon-rich (O-, C-rich), or even a combination of both. This is 
determined by the abundance of each element during the AGB phase. Spectroscopic observations did not reveal any 3.3/mi poly- 
aromatic hydrocarbons (PAHs, i.e. signature of C-rich chemistry) in the core of M2-9 ( |Smith & McLean"||2008| l, but other PAH 
features are present in the lobes (this work, Section 3.1). The nebula is O-rich with C/O < 0.5 ( Liu et al.|2001| l, thus the dust is 
expected to be O-rich. The model contained only amorphous silicate dust, since the exact percentage of crystalline silicates that is 
locked within the disk is unknown. 



We have chosen a certain position angle (PA.) of the gaseous nebula to fit the observed visibilities. According to Phillips & 
|Cuesta|fl999| l, the PA. of M2-9 is -2°. We support this estimation from the positions of the two ansae at ~52" from the core and 
finding the best fit at the same position angle (Table|3]l. 



4.2. Fitting results 

Following the success of fitting a disk model for the core of Mz3 and knowing that those two objects are similar, we chose parameter 
ranges for the disk in M2-9 similar to those of Mz3 (Table[3]). A more in-depth presentation on general MC3D fitting, can be found 
in |Sauter et aL] l |2009) . 

- The disk's inner rim radius was chosen by fitting the observed visibilities with Gaussian uniform disks, thus estimating their 
diameters in milliarcseconds, and converting those into physical units. 



- We adopted an initial value for the mass of warm dust derived by Smith & Gehrz ( 2005 1, but it was increased in order to recreate 



the silicate absorption features. Parameter a was fixed likewise to a certain value, where we would see silicates in absorption 
and at flux levels between those set by ISO and MIDI spectra. In this paper, we do not probe the large amount of cold dust that 
was observed by Smith & Gehrz (2005 ). 
- Low visibilities indicated a resolved structure, or else a very dense medium, so a scale height, ho higher than 20 AU, was chosen 
initially, with increments of 5AU, and a moderate flaring, j3, with steps of 0.01, in order to fit the observed visibilities (Fig. [7]). 

With the use of MC3D, the reconstruction of the disk's image was made. The models of its spectral energy distribution (SED) 
and visibility amplitudes have been compared to the MIDI data and the observed SED. After careful examination of all possible 
parameters, the model with the best fit, is one that resembles the best geometrical shape for a flared disk. The model's parameters 
are shown in Table [3] 

We were able to define a lower and an upper limit to constrain the size of the disk's structure from a theoretical basis and from 
our observational data. A projected baseline on the equatorial direction of M2-9 would probe the full size of the disk. In our case, 
this telescope configuration was not accomplished. 

The distance of 640pc suggested by Schwa rzet al.| ( |199"7] ) led us to many difficulties in fitting the observations. Even with scaling 
down the luminosity and the size of the disk or modifying the dust content, thus the optical depth, no improvement was found in 
fitting both the SED and the visibilities. Significant geometrical structures should have been detected in the dispersed visibilities 
for such a nearby source, yet the observed dispersed visibilities are very smooth. None of our models could fit the distance and 
luminosity assumed so far. Our best fits were found for a larger distance, namely ~1200pc (Table [3] Figures [6] [7] and [8]). Recent 
calculations^Jput M2-9 at a distance of 1300 ± 200pc (R. Corradi, private communication). 

The mass of dust, density and geometric parameters used in this model, represented the best fit for the opaque structure detected 
by MIDI in our line-of-sight (74°). A similar inclination of 73° was also found independently by Solf ( 2000| l for the polar lobes and 
by |Zweigle et al.| ( |1997| l for the CO torus. 

To fit the near- infrared flux, we have selected a cooler star (15,000 K) as the illuminating source than what was proposed 



by |Smith & Gehrz (2005 1 (30,000 K). From Fig.|6]one can notice that the shape of both MIDI and ISO spectra is the same; the only 
difference is the flux level. We preferred to fit the ISO spectrum, because our model has a larger aperture than that of MIDI. Our best 
fit for the MIDI visibilities was found by optimising the model's total x 2 , taking into account the restricted field-of-view of MIDI. 
The current fit has a total x 1 = 3.17 (Fig. [Til. 

We preferred to fit the visibilities rather than achieve a better fit for the SED. We could not fit our data by decreasing the outer 
radius, since these modifications increased the silicate absorption feature and lowered the visibilities, suggesting either a denser 
or a larger, resolved structure (without altering the shape of the spectrum that was due to colder dust). Expanding the outer radius 



7 Revision of the expansion data. 
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Parameters 


M2-9 


Mz3 


Teff (K) 


15,000 


35,000 


Luminosity (Lq) 


2,500 


10,000 


Distance (kpc) 


-1.2 


-1.4 


Inclination (°) 


74 ± 1 


74 ± 3 


best PA. (°) 


-2 ± 2 


5 ± 5 


Inner radius (AU) 


15 ± 1 


9± 1 


Outer radius (AU) 


900 


500 


a 


2.0 ±0.1 


2.4 ± 0.1 




0.9 ±0.1 


1.02 ±0.02 


frlOOAU (AU) 


36 ± 2 


17 + 2 


p gmm (gem" 3 ) 


2.7 


2.7 


Dust mass (Mq) 


1.5 x 10~ 5 ± 5 x 10" 6 


9 x 10~ 6 ± 2 x 10~ 6 



Table 4. Accretion luminosities for a range of WD masses and different companion mass outflows. 



outflow 






L acc (Lq) 






(kms- 1 ) 


0.6 M© 


0.8 M© 


1.0 Mq 


1.2 M© 


1.4 Mq 


5 


1,144 


2,760 


5,460 


8,734 


14,740 


10 


73 


174 


338 


566 


914 



imposed a new problem: more cold dust must be deposited in the disk, thus the predicted flux of the model exceeds actual measured 
levels at wavelengths longer than 12yum. For our best fit, the disk's optical depth is in the N band is tio^ — 4.6. The fit is sensitive 
to parameters a, ft and ho, as small changes in those values instantly altered the disk's opacity. The model in use does not include a 
puffed-up inner rim. 

A single illuminating source is assumed for the dusty disk; in this case the red giant. Thus, the contribution from the companion 
in both the heating and shaping processes is not included and cannot be reproduced. This drawback may partly explain the difficulty 
in fitting the spectral energy distribution. The binary interactions lead us to consider a different disk geometry. 

Due to the limited number of baselines, our selection of geometric parameters was restricted. We have found an acceptable disk 
model to fit the observational data, which implies a general agreement in terms of emission and spatial distribution of the bulk of 
the dust in the core of the system. Despite this, it is also true that the source is complex, probably engulfing a binary, and that the 
passive disk model is probably not a good description of the reality. This source clearly requires interferometric imaging either with 
aperture masking techniques ( Tuthill et al.|2006|l or with MATISSE/VLTI (Lopez et al.|2008 ). 



5. The mass of the components and an accretion disk 



Spectropho tometry of the core has shown an intrinsic amount of (high-excitation) iron emission lines ([Fell] and [Fe HI], [Allen & 
Swings 1972 1. This gives an electron temperature for the core -10,000 K and an electron density of 10 7 crrT 3 (similar to Mz3, Smith 



2003 ). As it is mentioned in that work, at such high core densities collisional de-excitation prevails and weak forbidden nitrogen 
lines in the core are not a surprise. Iron was found in the northern lobe as well (Fig. |3j. H ora & Latter| (1994) have also detected 
iron lines that have a high ionization potential, evidencing the presence of a particle beam or supersonic shock interface inside 



M2-9. Ha emission from the core is optically thick with broad wings due to scattering ( |Balick|[T989| |Arrieta & Torres-Peimbert 
2003). Polarisation measurements in the optical revealed an elliptical polarisation in the core of M2-9, which was consistent in every 
observed waveband (namely B, R and I). Aspin & McLean ( 1984) explain this as scattered continuum from the central ionizing 
source by an extended ionized dusty torus around the core (diameter < 4"). All this indicates that the illuminating source of M2-9 
must be a hot and compact star. 

A lifetime of a planetary nebula is -10,000 years and that of a white dwarf billions of years. Thus, material ejected by what is 
now the secondary, a white dwarf, is already dissolved in the interstellar medium. We expect that the material that composes M2-9 
at the moment, originates from what is now the primary, an evolved giant, either an AGB or a post-AGB star. 

According to Akashi & Soker] ( |2008[ l shaping of the ejecta of a tight-waist PN is possible by an accretion disk around the 
secondary. An accretion disk would be bright enough and illuminate part of the disk that has been detected by MIDI. We believe 
that the disk found in the core of M2-9, engulfs an evolved star (AGB or post-AGB) and a WD, namely the primary and secondary 
components of the binary. The first one is the main heating source for the dust at -2,500 Lq (Table [3]), while the latter obtains its 
luminosity from an accretion disk, although it is the main ionizing source for the nebula (lighthouse beam). 



The most probable scenario to explain the observables of Schwarz et al. ( 1997 1 is that the mass-losing star transfers material 
onto the secondary at a rate of M* = 10~ 6 MQyr -1 at supersonic speeds (5-10km s" 1 ) feeding the accretion disk of a white dwarf. A 
luminosity of -100 Lq was estimated for the ionizing source from radio observations ( Purton et al.|19 82). The following equation 
by Soker (2004 1 was used in order to calculate accretion luminosities for the white dwan: 



3 x 10" 



WD 



0.6M o 



) (lOkms 1 ) 
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Table 5. Dust sublimation radii in AU for two stellar luminosities 



T S ub [K] 


Rsub 


[AU] 




2500 L Q 


100 Lq 


2000 


0.97 


0.19 


1500 


1.74 


0.35 


1000 


3.90 


0.78 



(— n 

U00AU/ \ 



■ (3) 



10- 4 M o yr 



where M acc is the mass accretion rate, M* is the mass-loss rate of the primary, Mwd is the mass of the accreting white dwarf, a is 
the binary separation and v s is the wind speed. 

Table |4] displays a range of accretion luminosities (L a cc=GM acc MwD/^WD) for different secondary masses for the above- 
mentioned accretion rates and velocities, assuming a binary separation of 30 AU and a WD radius of 6,000 km. Our findings 
confirm a WD mass of 0.6-0.8 Mq, previously suggested by |Livio & Sok er (2001]) and in accordance to |Gesicki & Z ijlstra ( 2Q07] > . 



We can see that the accretion luminosities for a wind velocity ~10kms 1 aren't high enough to evaporate dust up to 15 AU 
(inner rim). The cool primary's luminosity (2,500 Lq) is sufficient though. 



6. Where is the companion relative to the dust source? 

It has long been believed that M2-9 engulfs a binary system. The light-house effect is an indirect, but rather convincing evidence of 
binarity. The assumption of circularity is based on the fact that no evide nt acceleration nor dec eleration of the lighthouse's angular 
velocity has been observed to date (R. Corradi, private communicationf*j |Livio & Soker| ( |200T| l have used a typical mass for a white 
dwarf of 0.6 Mq and 0.8 Mq for a post-AGB star, leading to an orbital separation of 27 AU, assuming a circular orbit. 

The VLTI observations of a compact dusty structure impose some constraints, which are not trivial for the interpretation of 
this system. However, these interferometric observations are also limited in the sense that they do not provide any closure phase, 
nor absolute astrometry as in radio interferometry. It is therefore not possible to locate precisely the position of the compact dusty 
structure. 



6.1. The inner dusty disk. 

The size of the dusty structure detected by MIDI by the baseline closer to the equatorial plane (M2-9_l; Table [T} is approximately 
33 mas (Fig. [5]). Extrapolating to a projected baseline at 90° (equatorial plane) we would expect a minimum size of 41 mas. From 
that we presume that the central cavity must be smaller than this size. Our modelling has showed that a diameter of approximately 
30±2 AU (or ~25 mas) fits the inner portion of the disk (Table|3]l. We did not find a better fit to the MIDI data than the one presented 
here by increasing the disk's inner rim. 

We have used a relation for dust sublimation radius (R su b in AU) by Tuthill et al. ( 2001| l: 



R 



sub 



1.1 



/ T sub \- 
U500K/ 



\ IOOOLq 



(4) 



where L is the stellar luminosity and r su b is the dust sublimation temperature. For a star of 2500 Lq (Table[5]) it is suggested that the 
inner rim of the disk must reside at a distance larger than 4 AU but within the limits defined by the MIDI observations. An 100 Lq 
heating source, such as the accretion disk in the case of M2-9, does not fit the observational data; neither did the 553 Lq suggested 
by|Schwarz et al.|(|1997|. 



We estimated orbital periods for the binary by using the observations of 18 years from Doyle et al. (2000 1. From those, it can be 
seen that the lighthouse beam covered less than 20% of the perimeter of the lobes within that time scale. We found that the period 
ranges from 90 to 120 years. 

A binary system should reside within the structure detected by MIDI. By applying simple Keplerian physics (assuming circular 
orbits), we have used several mass ratio^jfor the binary components to estimate the corresponding orbital diameters, for periods of 
90 and 120 years, in physical units and converted them to angular sizes for a distance of 1.25 kpc (Table[6|). 

As seen in Table [6] orbital diameters that fit within the extrapolated angular size of 41 mas range from 42-52 AU. This is larger 
than our best fit model, where the size of the disk's inner cavity is 30 AU. We should draw attention though to the fact that the 
disk's inner rim is centered at the heating source, which in turn is not positioned on the binary's centre-of-mass. The model in use 
cannot reproduce a partially illuminated disk or a disk with an off-centre illuminating source. Thus, we believe that the companion 
is truncating the disk, since its heating capabilities are minimal. 

We can only surmise that high-angular resolution radio observations could establish the astrometric position of the companion 
(WD/accretion disk). 



Bearing in mind though, that intensive monitoring of the nebula occurred only in the last 20 years (Doyle et al. 
9 Within the Chandrasekhar limit for the WD 



2000) 
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Table 6. Binary orbital diameters in AU (also converted in mas for a distance = 1.25 kpc), at two different periods, 120 and 90 years 



Mi (Mq) 


M 2 (Mq) 


P = 


120 yrs 


P = 


90 yrs 






d[AU] 


d [mas] 


d[AU] 


d [mas] 


0.6 


0.6 


52 


41.6 


42 


33.6 


0.8 


1.0 


60 


48.0 


50 


40.0 


0.6 


1.4 


62 


49.6 


50 


40.0 


0.8 


1.4 


64 


51.2 


52 


41.6 


1.0 


1.2 


64 


51.2 


52 


41.6 


0.6 


1.0 


57 


45.7 


47 


37.7 


0.6 


1.2 


60 


48.0 


49 


39.2 


1.0 


1.0 


62 


49.2 


51 


40.6 



6.2. The exterior torus 



The obtained NACO images are fully dominated by the central source and lack dynamical range in its close vicinity (Fig. [TJ. At 
~3" from the core |Zweigle et aL] ( |1997 1 found that the CO emission is coming from a large ring, whose center seemed offset 
by 0.5+0.3" from the compact radio source (< 0.1" and unresolved at 1.3 cm, Bignell p983| ). More recent Plateau de Bure data 
(Castro-Carrizo et al., in preparation) with 0.3" angular resolution confirm the presence of an offset, and detect additional CO 
emission coming from regions at a distance of ~1" from the center. 

Molecular hydrogen spectroscopy has yielded two different velocity components at distances <0.5" latitudinally fro m the central 
illuminating source: 10.9km s _1 (blue-shifted) and 123.7km s -1 (red-shifted). This points to a H2 disk-like structure (Smith et al. 
[20051 . 

Spiral structures (as expected in a symbiotic type system) can be formed by binary interactions during the ejection of the 
primary's CO envelope (Edg ar et al.|2008| , yet this is not clear at this point for the case of M2-9. If there are spirals, they are seen 
at a high inclination (pole-on). Our source is edge-on. The differential phases of the disk detected by VLTI are small, lying within 
a range of +10°, i.e. well below some strong signature from dusty rings reported in the literature ( Deroo et al.|2007a Ohnaka et al. 
[2008] l. 



7. Comparing disks: M2-9 and Menzel 3 

These two bipolar nebulae are spectroscopic twins and both exhibit complex gaseous shell structures. Both are extended at similar 
sizes. Their main parts, including the lobes, are 43" and 38" at length and ~15" at width, for M2-9 and Mz3 respectively. The ansae 
of M2-9 are 115" apart, while the farther regions of Mz3 are at about 130". 

Distances to both objects are not well determined and the estimates given stand for the best fits to the observational data: M2-9 at 
1 .2kpc, Mz3 at 1 .4kpc. Indirect evidence for binarity in the case of M2-9 are the lighthouse knots in the nebula and the fast outflows 



observed in Ha (as recently studied in the case of the Red Rectangle, Witt et al. (2009 1), and for Mz3, the X-ray jet (Kastner et al. 
|2003| l and also fast outflow seen in Ha. 

Whereas Mz3 was unresolved in all infrared filters, some significant extensions of the nebular core at larger diameters in both 
equatorial and polar directions were found for M2-9, as already suspected in some other observations. Our observations in the 
mid-infrared have revealed that they both have very compact dusty disks in their cores. 

The disks have been resolved with high-angular resolution techniques and their composition is similar (amorphous silicates). In 
addition, the disks are in the equatorial planes of the planetary nebulae (that is, perpendicular to the bipolar ejecta) and at similar 
inclination from the rotation axis (~74°). 

The disk inside M2-9 is larger in terms of size (Fig. Ml and hot dust mass compared to Mz3: 1.5 x 10~ 5 M0 and 9 x 10~ 6 M0, 
respectively. In both cases, the total amount of materialis 0.5-1 Mq, which points to progenitors of several Solar masses. This 



is in line with the high nitrogen content found by Smith (2003) for Mz3. In addition, the disk of M2-9 exhibits some degree of 
crystallinity (Fig. |4), which may indicate that it is older, or that the physical conditions for annealing prevail. There are no such 
features in Mz3 (Fig. [9]). 

In the work of S mith| ( |2003 1 > it is shown that the optical to infrared SED of Mz3 required two components; a hot star and a cool 
star. M2-9 has a very similar SED in the optical to that of Mz3, thus it is probable that it has similar components. 

Assuming a circular orbit might be an oversimplification. Binaries containing a WD and another evolved star, may have eccentric 
orbits. When main sequence stars will ascend to the red giant branch and then to the AGB, it should be expected that despite strong 
mass-loss the eccentricity will persist, leading to variations in the building up of an accretion disk and the generation of jets. This 
may lead to short and strong events of large mass ejection, related to the building of large scale structures seen in the nebula, and in 
the equatorial plane as well. 



8. Conclusions 

We present the discovery of a complex dusty disk inside the elongated bipolar planetary nebula M2-9, with the use of infrared 
interferometry, for which we have derived geometrical constraints and chemical composition. Our findings show that the disk is 
much smaller than the extended dusty structure around the nebula's core with an inner rim of 15AU at 1 .2kpc and that it is primarily 
composed of silicates. This suggests that the dusty material is derived from the AGB envelope of the primary and has survived 
its evolution onwards. The disk is also aligned with the equatorial axis of the nebula. It is probable that the dust is heated by the 
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primary, evolving star (O.6-1.4M0), but it is also partially heated and truncated by the accretion disk of the secondary, a white dwarf 
(0.6-1 M ). 

There are similarities between the disk found in M2-9 and that of its spectroscopic twin, Mz3. Although the sample is small, 
many bipolar planetary nebulae and proto-planetary nebulae (e.g. M2-9, Mz3, OH231.8+4.2, He2-1 13, CPD-56°8032) have dusty 
structures (disks/tori/spirals) around their cores. Some of them are known to engulf binary systems. Binarity in the cores and dusty 
structures around them may explain the shaping of the ejecta into bipolar shapes. Thus, an extended study of those bipolar proto- 
planetary and planetary nebulae would be advised to establish whether this is a global tendency. 

M2-9 is yet another example of an object with silicate chemistry within the disk, suggesting that the material was ejected during 
the oxygen-rich phase of the primary star. Nevertheless, there are PAHs in the lobes of the nebula, at a much lower percentage, most 
probably formed after the dissociation of carbon monoxide. 
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Fig. 3. High-resolution spectra of the northern lobe of M2-9 by IRS/SPITZER. Upper curve corresponds to the Short-High slit and 
lower curve to the Long-High (for slit orientation, see Pig. |2j. The inset in the upper plot is a close-up of the PAH emission features. 
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Fig. 4. Top: Normalized visibility vs. wavelength for M2-9 revealing two crystalline silicate absorption bands at 9.8 and 11.4yum; 
they are identified as forsterite. For the sake of clarity only the nearest to the equatorial plane measurement (M2-9_l: 40.2m, 107.3°) 
is shown here. Bottom: Mass absorption coefficients for forsterite (solid) and enstatite (dash) taken by Koike et al. ( 2000 2006) 1 
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Fig. 5. Gaussian size distribution of the disk (±5%) for each baseline (Table [ljl in N band implementing the dusty structure's 
evolution in the mid-infrared. 
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A (yum) 



Fig. 6. At 1.2kpc: MIDI spectrum (blue) compared with the ISO spectrum (black - solid: 1 4" x20" aperture, dashes: 



14"x27"aperture), ISOCAM data (green crosses), AKARI data (open circles), the TIMMI2 photometric measurements of Smith & 
Gehrz (2005 ) with an aperture of 4" (cyan triangles), IRAS measurements (circles) and Cohen & Barlow ( 1974| > photometry (purple 



stars, 1 1" beam). The NACO L' and M' band are also included (magenta stars, aperture 0.4")- The red line is the best model with 
the full aperture (~1")- All spectra have been dereddened using the Savage & Mathis ( 1979| law. 
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Fig. 7. MIDI visibilities, compared with the best model (dashed-dotted lines, ^ 2 =3. 17, see Table[3J». The visibility curves are for the 
baselines M2_9-l to M2_9-4 from top to bottom and left to right (Table[TJ. Individual^- 2 are 2.61, 5.33, 3.05 and 1.69 for each curve 
respectively. 
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Fig. 8. Flux distribution of the model at 8, 10 and 13 /zm. Images on the top belong to M2-9 and the ones on the bottom to Mz3. A 
small amount of light emerges at 8 yum while the flux comes from inner disk parts at 13 /urn. 
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Fig. 9. Normalized visibilities vs. wavelength for Mz3 ( |chesneau et al. 2007 1. The prominent feature belongs to the bright 12.8//m 
[Nell] nebular emission line. The noise increase at 9.1pm corresponds to the narrow ozone feature. Observed baselines are indicated 
within the legend. 



